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Abstract 

We propose a new non-perturbative approach to hadronic interactions at high ener- 
gies and small momentum transfer, which is based on the scale anomaly of QCD and 
emphasizes the role of semi-classical vacuum fields. We find that the hadron scattering 
amplitudes exhibit Regge behavior and evaluate the intercept o(0) of the corresponding 
trajectory. Both the intercept and the scale for the slope of the trajectory appear to be 
determined by the energy density of non-perturbative QCD vacuum (the gluon conden- 
sate). Numerically, we find A = a(0) — 1 = 0.08-^0.1, consistent with the values ascribed 
phenomenologically to the "soft" Pomeron. For arbitrary numbers of colors N c and flavors 
Nf, A is found to be proportional to (Nf/N c ) 2 ; however, in the large N c (Nf fixed) limit, 
A ~ N°. 



Understanding the behavior of QCD at high energies and small momentum transfer is still 
a challenging and unsolved problem. In the framework of perturbation theory, a systematic 
approach was developed by Balitsky, Fadin, Kuraev and Lipatov |[|, who demonstrated that 
the "leading log" terms in the scattering amplitude of type (g 2 In s) n (where g is the strong 
coupling) can be re-summed, giving rise to the so-called "hard" Pomeron. Diagramatically, 
BFKL equation describes the t— channel exchange of "gluonic ladder" ( see Fig. la ) - a concept 
familiar from the old-fashioned multi-peripheral model @. 
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At small momentum transfer, QCD perturbation theory is in general inapplicable, but one 
may choose to consider the scattering processes where the parton virtualities at the ends of the 
ladder are fixed to be large (for example, the scattering of two heavy quarkonium states @). 
However, even then the partons can still "diffuse" to small values of transverse momenta toward 
the center of the ladder (diffusion in the log of transverse momenta |I|]), and at sufficiently high 
energies the perturbative approach inevitably breaks down [|]]. This argument was formulated 
rigorously by A.H. Mueller ||, who showed that the operator product expansion (which provides 
the basis of the perturbative approach) breaks down at high energies. Another serious problem 
of the perturbative treatment has been made apparent by recent vigorous calculations of the 
next-to-leading corrections to the BFKL equation ||. The NLO corrections appeared to be 
large, and drove the intercept of "BFKL Pomeron" significantly below the range of values 
suggested by phenomenology. 

Perturbative expansion of the scattering amplitude is possible only in the presence of a 
sufficiently large scale. As was mentioned above, at very high energies the external scale, which 
determines the parton virtualities at the ends of the ladder, becomes progressively unimportant, 
and the perturbative expansion loses justification. Therefore it looks plausible that the Pomeron 
is a genuinely non-perturbative phenomenon 0, ||. At present, non-perturbative phenomena 
can only be treated theoretically if they stem from relatively short distances, which requires the 
presence of a large scale. The main idea exploited in this letter is that such a scale exists in the 
QCD vacuum as a consequence of scale anomaly, and is related to the density of vacuum gluon 
fields due to semi-classical fluctuations; numerically, Mq ~ 4-=-6 GeV 2 (see below). Because of 
the presence of this large scale, the perturbative expansion still makes sense; we are able also 
to evaluate explicitly the leading non-perturbative contribution due to the scale anomaly. 

There are two facts that support the feasibility of such an approach. First, the success of 
QCD sum rules is based on the use of a few first terms in the operator product expansion, which 
can be justified only if a sufficiently large scale associated with the vacuum structure exists ||. 
Second, the non-perturbative amplitude of low-energy dipole-dipole scattering was evaluated 
and found to be determined by the vacuum energy density, arising from the semi-classical 
fluctuations of gluon fields |[L0|| . This latter example is encouraging, since the multi-peripheral 
model 0] relates the amplitude of high-energy scattering to the low-energy interactions of 
partons. 

Basing on these ideas, we propose an extension of the BFKL program to the non-perturbative 
domain; the key ingredient of our approach is the breakdown of scale invariance in QCD, re- 
flected in scale anomaly. The concept of scale anomaly is rather general and was formulated 



long time ago [ ] T|,[12 ],[13|; let us briefly recall its application to QCD. In the chiral limit of 
massless quarks, the Lagrangean of QCD is scale invariant on the classical, tree level. This 
invariance is however broken by renormalization, which introduces a dimensionful scale once 
the interactions are switched on. This "dimensional transmutation" phenomenon is fundamen- 
tal for the understanding of scale dependence of the strong coupling constant |14|], which is 
the basis of all applications of perturbative QCD. On the formal level, the breakdown of scale 
invariance in the theory is reflected by the non-conservation of scale current, and thus in the 
non-zero trace of the energy-momentum tensor 9^ u [15]. Scale anomaly leads to a set of pow- 
erful low-energy theorems for the correlation functions of gluon currents in the scalar channel 
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B- 

The starting point of the approach that we propose in this letter is the following: among the 
higher order, 0(a 2 s ) (as = g 2 /47t) , corrections to the BFKL kernel we isolate a particular class 
of diagrams which include the propagation of two gluons in the scalar color singlet channel 
J PC = ++ (see Fig. |l]-b). We will show that, as a consequence of scale anomaly, these, 
apparently 0(a|), contributions become the dominant ones, 0(a s ). 




Trvmnmrr 

colliding hadrons gluons pions 



Figure 1: Multi-peripheral ( ladder ) diagrams contributing to the leading-order BFKL (a) and 
"soft" (b) and (c) Pomeron structure 

Indeed, let us consider the contribution of Fig. Q-b, which is one of the numerous corrections 
of the next-to-leading order to the BFKL Pomeron. In perturbation theory, such corrections 
are oc a s ; however, we note that if the two produced gluons in Fig.|I]-b are in the scalar and 
colorless state, the vertex of their production, generated the four-gluon coupling in the QCD 
Lagrangean, is ~ a s F^ u a F® u . We observe that this vertex is therefore proportional to the trace 
of the QCD energy-momentum tensor ( 9fi) in the chiral limit of massless quarks: 



n/J, _ 0{g) paa/3 pa ^ ^9 paa/3 pa . f-t \ 

"~ 2g a(3 ~ 32^ Q/3 ' 



note that as a consequence of decoupling theorem |L6| the (3 function in Eq. (0) does not contain 
the contribution of heavy quarks (i.e. b = §(lliV — 2Nf) = 9). 

The entire contribution of Fig.|TJ-b therefore appears proportional to the correlator of the 
QCD energy-momentum tensor. Let us now consider the spectral representation for this cor- 
relator: 

IT(g 2 ) = i Jd*x (0|T {0*(xW(0)} |0> = / da 2 -^^— , (2) 
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with the spectral density defined by 



p e (k 2 ) = Y,(^m Pn -k)\(n\6Z\0)f 



(3) 



where the phase-space integral is understood. In lowest-order perturbation theory, the spectral 
density (|J) is given by the contribution of two-gluon states; the calculation for SU(N) color 
gives 



bg 2 V N 2 



32vr 2 



471-2 



(4) 



However, at small invariant masses, perturbation theory inevitably breaks down; an impor- 
tant theorem || for this correlator states that as a consequence of broken scale invariance of 
QCD, 

n(O) = -4<O|0£(O)|O). (5) 

Since this theorem, as will become clear soon, is a corner-stone of our approach, let us 
briefly recall its proof ||. It is based on the fact that the expectation value of any operator O 
of canonical dimension d (d = 4 for 9fi) can be written down as 



(O) 



M exp 



57T 



bgl 



(6) 



where go = g(Mo), and Mq is the renormalization scale. On the other hand, the dependence 
of QCD Lagrangean on the coupling is (— 1/4<7q) F aa ^F^g, where F = gF is the rescaled gluon 
field. By writing down the expectation value of the operator O in the form of the functional 
integral, and by differentiating this expression with respect to 1/4^q, one can therefore generate 
correlation functions of the operators O and F 2 . Differentiating once, one gets 



dx (T{0(x) F 2 (0 



d 



■(O). 



(7) 



Combining (|7D and (j^), and choosing 0(x) = 0^(x), furnishes the proof of the theorem (|5|). 

Note that the r.h.s. of Eq. (||) is divergent even in perturbation theory, and should therefore 
be regularized by subtracting the perturbative part. The vacuum expectation value of the 9£ 
operator then measures the energy density of non-perturbative fluctuations in QCD vacuum, 
and the low-energy theorem (||) implies a sum rule for the spectral density: 



da 2 



a- 



[P 



phys (a 2 ) 



pgV)] = -4 (0|^(0)|0> = -16 e vac ^0, 



(8) 



where the estimate for the vacuum energy density extracted from the sum rule analysis gives 
e vac ~ —(0.24 GeV) 4 |17| . In addition, another sum rule 



18, 17 



da 2 p^ s (a 2 ) 



da 



2 ^(a 2 ) 



(9) 
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is implied by the quark-hadron duality. Since the physical spectral density, Pg hyh , should 
approach the perturbative one, pjf, at high a 2 , the integral in Eq. (||) is convergent. 

According to (|l|) and (Q), the l.h.s. of Eq.(|[) is apparently 0(g 4 ); however it is easy to see 
that this is not so by looking at the r.h.s. of this equation, which is renormalization group 
invariant, and does not depend on the coupling constant. This means that the l.h.s. must 
also be 0(g°). Let us illustrate this formal argument by considering the spectral density (|3|) at 
small invariant mass fl9| . Small invariant masses imply small relative momenta for the produced 



particles, and at small momenta an accurate description of QCD is given by an effective chiral 
Lagrangean 

£ = & tr d^Ptf + i mlfl tr (u + , (10) 

where U = exp (2iir/f n ), it = n a T a and T a are the £77(2) generators normalized by tr T a T b = 
^5 ab . The trace of the energy-momentum tensor for this Lagrangian is (see, e.g., [[TTJ) 

f 2 

e» = -2 ^ tr d^U^rf - m 2 Jl ^ (u + U ] ) . (11) 

Expanding this expression ( |TTj) in powers of the pion field, one obtains, to the lowest order, 

0£ = -d^dV + 2m 2 7r a vr a + • • • , (12) 



and this leads to an elegant result |19| in the chiral limit of vanishing pion mass: 

<7r+7r-|^|0> = g 2 . (13) 

This result for the coupling of the operator 8£ to two pions can be immediately generalized 
for any (even) number of pions using Eq. (|TTJ) . The expression (|13"1) is manifestly ~ 0(g°), 
and shows that the spectral density of the scalar gluon operator ~ g 2 F 2 is independent of the 
coupling constant g as a consequence of scale anomaly. While we have used an effective chiral 
Lagrangean to illustrate how the dependence on the coupling constant gets "eaten" by the 
scale anomaly, this phenomenon is very general and does not depend on the specific model for 
the spectral density. One way of understanding the disappearance of the coupling constant in 
the spectral density of the g 2 F 2 operator is to assume that the non-perturbative QCD vacuum 
is dominated by the semi-classical fluctuations of the gluon field. Since the strength of the 
classical gluon field is inversely proportional to the coupling, F ~ l/<?, the quark zero modes, 
and the spectral density of their pionic excitations, appear independent of the coupling constant. 

Armed with this knowledge, we are ready to see that the contribution to the next-to- 
leading order BFKL kernel that describes the production of two gluons in the color singlet, 
scalar state, which is formally ~ 0{g i ) 1 as a consequence of scale anomaly can become the 
leading one, ~ O(g ). (Of course, the perturbative part of this contribution is still ~ 0(g 4 ) 
and has been taken into account in the next-to- leading order BFKL Pomeron). Therefore, we 
want to build a multi-peripheral model for the "soft" Pomeron in which hadrons are produced 
(mostly two pions, see Refs.|TD|, |2"0]] ) due to exchange of two gluons in the t-channel (see 
Fig.|l]-c ). The only dimensional scale in this approach appears in fP hys and can be estimated 
directly from sum rules of Eq. (H). It turns out that the characteristic mass (Mq) in Eq. (§) 
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is rather large [Kj , M5 2 « 4 GeV 2 (the original analysis of |J , [^TJ yielded even somewhat 
bigger value Mq rs 6 GeV 2 ). This is the largest scale which exists in non-perturbative QCD 



22]. In the framework of the instanton approach, the large magnitude of Mq was shown to 



be a consequence of strong color field inside the instanton [22 



This value determines the scale of all dimensional parameters of the Pomeron trajectory as 
well as the typical transverse momentum of produced particle in the Pomeron. It is interesting 
to notice that the experimental value for the slope a' P (0) of the Pomeron trajectory (in the 
standard notation, a P (t) = 1 + A + a' P (0) t) a' P (0) = 0.25 GeV~ 2 [g3|] is very close to 1/M 2 . 

We start our calculation with diagrams of Fig.|l]-b in the leading log s approximation of 
pQCD where we sum only contributions of the order of (as lns) n . In this approximation the 



propagators of the t-channel gluons can be written in a simple form [24 



^tiAQi) - — 2~ - ~2~ x 

Hi Qi,± 



CXif3iS 



(14) 



where we use the Sudakov decomposition for momenta qi along the momenta of colliding par- 
ticles ( pi and P2 in Fig.|l]-b ), namely, 



(15) 



Eq. ([14]) corresponds to Weizsacker- Williams approximation for the gluon field of a fast-moving 
hadron. 

The ladder diagram of Fig.|l|-b for emission of n-pairs is equal to 



a n (Q 2 ) 



a 



2 f 1 f^ 1 s daid^d 2 q i!± 
sJliFv 11 o/-o-\a 5(aif3 i+1 s 

i=l 



2(2tt) 



ML) *i 



(16) 



r^ <j ^ i+1 r fc , ijI/j+:L G^ Wj ^ <+1 (^ )Gi/ u v i+1 {{Q Qi) ) r 



where cr(Q 2 = 0) is the total cross section of n-pairs production (cr n ) and $j is the phase space 
factor for two identical particles with total mass Mj which is equal to 



d 3 ^ 



(2vr) 3 2o;i2a;; 



S(ui + u> 2 — Mi 



327T 2 



(17) 



In Eq. (0) r MiiA t i+1 is the vertex of gluon pair production. It is easy to calculate that it is equal 
to 

Qi,±,fJ,i Qi,±,lH+i ^fH,IH+l = 3 %x • $+1,-1- , (18) 

after projecting on the colorless state with J PC = ++ . 

Using the kinematic relation a, s = Mf + kf x , where k^± = — $+i,± and per- 
forming integration over explicitly, we can rewrite Eq. ([16]) in the simple form 



a, 



(In s 



n i=n+l 



n 



4 • 9 • a 2 s (Q± - q it± ) ■ (Q± - q i+ljl _) 



dM} 



32tt 5 



qh (Q - q;" 2 



(m 2 + kl 



(19) 
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For forward scattering Q\ = and Eq. (|T9| ) leads to power-like behavior of the total cross 
section: 

oo 

\^ BORN A / on \ 

&tot = O n = O S , (20) 

n=0 

where 

al 18 /■ d*; 2 : rfM 2 



32tt 2 7 (M 2 + klf V ; 



and the cross section due to two gluon exchange 



a BORN = a 2 



(22) 



Eq. ( pi] ) can be easily rewritten through the perturbative spectral density p p e QCD that was 
evaluated above (see Eq. (f|)): 

where a 5 (M 2 ) = 4vr/(61n(M 2 /A 2 )). 

Our main idea is to separate non-perturbative and perturbative contributions to the spectral 
density of the scalar gluon operator by replacing p p e QCD {M 2 ) with (p P 1 hys (M 2 ) — pp® 00 (M 2 )) + 
Pq (M 2 ). The purely perturbative contribution is of the order of 0(g 4 ), and has been 
evaluated before ||. For the non-perturbative contribution, in which we are interested here, 
we have 

A = T x (t) 2 x 3^ f % (^'< M ' 2 > " ^ CD < M2 > ) ' < 24 > 

To estimate the integral in Eq. ( p4|) we will use the chiral approach to p p e hys described above, 
namely, 

pf ys {M 2 ) = ^- 2 M\ (25) 

62,71 

which corresponds to diagram of Fig. |l]-c Q. 

It is instructive to establish a qualitative relation between the matching parameter Mq and 
the energy density of QCD vacuum using the spectral density ( [25]) and the sum rule (§). Since 
perturbative spectral density pg® CD (M 2 ) at moderate M is much smaller than pg hys (M 2 ), Eqs 



and (pf) lead to the following approximate relation [21 



Ml c 327rj^_j , (26) 

which shows that the matching scale Mq is directly determined by the energy density of the 
vacuum. Since 

^vac ~ iVj, the magnitude of iVfo is proportional to N 2 jN 2 . 

1 Of course, the chiral approach cannot be extrapolated up to Mq ~ 2 GeV; at large M the spectral density 
(E9) should be corrected by a phenomenological form-factor expressed in terms of experimental irir phase shifts 
jbj; however, since our sensitivity to the large M region is only logarithmic (see (0)), in this paper we will use 
the simplified ansatz (|25|). 
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Collecting all numerical factors and substituting b = 9, we obtain 



A=lln4. (27) 
48 Ami 

Let us discuss the result of these simple calculations: 

1. Eq. (|20|) and Eq. (|27|) say that our approach leads to the exchange of the Pomeron with 
the intercept a(0) = 1 + A > 1. To the best of our knowledge, this is the only approach 
in the framework of QCD which leads to a "soft" Pomeron; 

2. Eq. (^) is a consequence of a direct generalization of the BFKL approach to the non- 
perturbative domain. It should be recalled that the only theory where the Pomeron 
naturally appears is the BFKL Pomeron in 2 + 1 dimensional QCD [p5 |; 



Eq. fl2"T|) gives A = 0.082 for Mq = 4 GeV 2 [|10j in good agreement with the phenomeno- 



logical intercept of the "soft" Pomeron, A = 0.08 [23]; it should be noted however that 



the precise value of the matching scale Mq as extracted from the low-energy theorem 

(HD depends somewhat on detailed form of the spectral density, and can vary within the 
2 _ a • a n„\r2 



range of M 2 = 4 -h 6 GeV 2 @, (ItJ. Fortunately, the dependence of Eq. (|27|) on M is 
only logarithmic, and varying it in this range leads to 

A = 0.08 -=-0.1. (28) 

As we have already stressed, in our approach the only dimensionful parameter is M^; 
its large value implies the dominance of rather short distances in the "soft" Pomeron 
structure. This fact is in agreement with a number of experimental and phenomenological 
observations: 

• The value of the slope for the "soft" Pomeron trajectory a' P (0) = 0.25 GeV~ 2 <C 
a' R (0) = 1 GeV~ 2 , where a' R is the slope of the Reggeon trajectory; 

• The experimental slope of the diffraction production of the hadron system with large 
mass is approximately two times smaller the slope for the elastic scattering. It means 
that the proper size of the triple Pomeron vertex is rather small. For our Pomeron 
it should be on the order of 1/M 2 « 0.25 GeV~ 2 < B el = 10 GeV" 2 ; 



The HERA data [26 on diffractive J/ty production in DIS show that the t - slope for 



elastic diffractive dissociation (7* +p — >• J/\P +p) is larger than the t - slope for the 
inelastic one (7* + p — > J/^f + X, where A is a high-mass hadronic system). This 
shows the existence of two different scales in the proton, one of which is determined 
by its size, and another one by the correlation length ~ 1/Mq of the gluon field 
inside. 
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5. Non-trivial azimuthal dependence observed recently in diffractive production of scalar 

if one adopts the idea that the effective coupling of the 



mesons [27] can be explained 



Pomeron to mesons is dictated by scale anomaly. 

The disappearance of the dependence on the coupling constant, which is the central point 
of our approach, may seem puzzling. However let us mention again that this result can be 
easily understood if we recall the interpretation of the non-zero v.e.v. of the trace of the 
energy-momentum tensor as being due to the semi-classical fluctuations of gluon fields. Since 
the strength of the classical gluon field is F 2 ~ l/cus, quark zero modes, and their pionic 
excitations, appear independent of the coupling, ~ 0(a^). We therefore envision the Pomeron 
as a t— channel exchange of two gluons, which scatter off semi-classical fluctuations of vacuum 
gluon fields; this scattering is accompanied by the excitation of quark zero modes in the vacuum, 
resulting in the production of pions. Amazingly similar picture of the "soft" Pomeron was 
anticipated by Bjorken ]3^]. These effects also manifest themselves in the low-energy scattering 
of heavy quarkonia; the magnitude of the scattering amplitude was found |TIJ to be determined 
by the energy density of the non-perturbative QCD vacuum. 

Let us discuss the dependence of our result on the numbers of colors, N c , and flavors, Nf. 



Two limits are of theoretical interest: i) N c —> oo, Nf, g N c fixed; ii) N c - 
fixed. The case i) corresponds to the large N c limit proposed by 't Hooft 
basis of "topological expansion" suggested by Veneziano M. 



oo, Nf/N c , g 2 N c 
while ii) is the 



Since the number of Goldstone bosons contributing to the non-perturbative spectral density 
for spontaneously broken SU^Nf) x SUii(Nf) is equal to Nj — 1, it is evident from Eq. (^) 

- 2Nj) in the denumerator) . A 
| Therefore, our approach may 



that A ~ N 2 f /N 2 (note that Eq. (g3j) contains b = 1/3(1UV C - 
simple graphic illustration of this dependence is given in Fig| 




oc N 



Goldstone bosons 



Fig. 2-a 




TTYYYmiT 



For large number of colors N c 

Fig.2-b 



oc N, 



Figure 2: A simple illustration of the appearance of the N 2 /N 2 factor in the Pomeron intercept. 

be considered as a realization of general ideas, proposed long time ago by Veneziano ||, that 
the "soft" Pomeron should be found keeping N 2 /N 2 fixed. 

Let us discuss now the large N c limit i), which corresponds to pure gluodynamics. Naively, 
since A was found proportional to N 2 /N 2 , one may conclude that in this limit A vanishes, and 
the cross section does not grow with energy. This conclusion is, however, immature. Indeed, 
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the physical spectrum in the scalar channel in gluodynamics contains a scalar glueball Q; its 
couplings to mesons are suppressed by l/N c (see, e.g., 0), so it should be very narrow. 
Therefore, the spectral density Eq. ( p5|) should be replaced by 

pf ys (M 2 ) = RM R 5(M 2 - M 2 R )+pert. contribution (29) 

where Mr is the scalar glueball mass, and R is its residue; the factor M R is introduced to make 
R dimensionless. Using Eq. (|29|) in the sum rule (H), we get a simple relation 

«=16^- (30) 



With Eq. (H) and Eq. ©, Eq. (||) becomes 

288^ M. 

since M R ~ iV?, e vac ~ AT?, and b ~ N c , Eq. (|31~D is well-defined in the large N c limit. 

One could try to estimate the value of A given by Eq. (|3lf) for Af c = 3; this requires the 



knowledge of the mass of the scalar glueball in pure gluodynamics. Recent lattice result |]3T 
gives Mr ~ 1.65 GeV; assuming that the main contribution to the energy density of the vacuum 
is due to gluons and using, as before, e vac ~ (0.24 GeV) 4 , b = llN c /3 = 11, we get the value 



^gluodynamics 



0.01, (32) 

which is significantly smaller than our result fl28|) for the world with light quarks. This indicates 
that the presence of light quarks in the theory leads to a much faster growth of the cross section 
with energy. 

The key question is whether one can prove the theoretical self-consistency of our approach. 
Indeed, our classification of the contributions to the scattering amplitude is still based on 
the expansion in powers of as, in which we have isolated the term ~ 0(a° s ) emerging as a 
consequence of scale anomaly. This term is the leading one only if the coupling constant as 
is sufficiently small. The magnitude of the coupling depends on the renormalization scale M . 
Since this dimensionful scale extracted from the sum rule analysis is large, M = 44-6 GeV 2 , the 
coupling constant indeed appears to be small, as{M^) -C 1. This fact insures that perturbative 
corrections to the kernel are smaller than the leading, ~ 0(a s ), term, and should be taken 
into account in the framework of conventional BFKL approach. Our approach yields a natural 
rapidity scale (Y = Ih(Mq / '4m 2 ) ~ 2 4- 3) for BFKL kernel above which the perturbative 
approach can be applied. It is interesting to note that for Y > Y the next-to-leading order 



corrections are well under control [ 33| ; however the interplay between "soft" and "hard" physics 
still has to be understood. 

We therefore believe that our proposal can lead to a systematic theoretical approach to 
the Pomeron in QCD. Of course, one cannot exclude a priori a different view of the Pomeron 



2 We thank S. Nussinov and E. Shuryak for stressing the role of the scalar glueball for our approach in the 
case of pure gluodynamics. 
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structure as coming from large distances, R 3> 1/Mq. However, experimental data (see point 4 
of our discussion above) suggest that the phenomenological Pomeron indeed originates at small 



distances R ~ yap(O) ~ 1/M . 

Let us discuss the relation of our approach to other existing approaches to soft scattering. 
Very similar ideas of the dominance of semi-classical vacuum gluon fields were developed in 



Refs. ||34fl . Our approach is complementary to these ideas, giving a natural explanation of 
the energy behavior of the soft scattering amplitude, which previously had to be taken phe- 
nomenologicallyQ. Let is note also that the correlation length of gluon fields, which was taken 



from the lattice QCD calculations in Refs [5^], in our approach appears only implicitly and is 
determined from the analysis of low-energy theorems. 

The dominance of classical gluon field configurations in high-energy collisions is the key 



idea of the approach proposed by McLerran and Venugopalan |36|] and developed in Refs. [[37] 



In this approach, the role of dimensionful parameter is played by the density of color charges 
in the transverse plane, rather than by the vacuum energy density. In our opinion, this is a 
plausible assumption at very high energies and/or for sufficiently heavy nuclei for low partial 
amplitudes (central region in the impact parameter plane). Since we focus our attention on 
the behavior of the total cross section, which is determined by large distances in the impact 
parameter space, and therefore small density of the color charge, the relevance of the scale Mq 
associated with the vacuum field strength should not be surprising. We feel that the approach 
of p6[ , [p^ j can describe the inclusive cross section, while ours is suited for the description of 



the total cross section. Indeed, the multiplicity associated with our multi-peripheral ladder 
~ 2 a(0) In s is rather small compared to the expectations of |36| , |37| . It would be extremely 



interesting to understand better the relationship between the two approaches. 
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